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(Dated: February 28, 2008) 

Abstract 

Using a 281 pb _1 data sample collected at the ^(3770) resonance with the CLEO-c detector at 
the Cornell Electron Storage Ring, we report the first observation of D + — > r]e + u e . We also set 
upper limits for D + — > r]'e + u e and D + — > (pe + u e that are about two orders of magnitude more 
restrictive than those obtained by previous experiments. 
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The quark mixing parameters are fundamental constants of the Standard Model (SM) 
of particle physics. They determine the nine weak-current quark coupling elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix Charm semileptonic decays have been 
studied in considerable detail because they provide direct measurements of the magnitudes 
of the CKM elements V c d and V cs , and a stringent test of theoretical predictions of strong 
interaction effects in the decay amplitude. In order to gain a complete understanding of 
charm semileptonic decays it is important to study as many exclusive modes as possible. 
Several rare modes have yet to be observed. 

The most precise measurements to date of absolute branching fractions of Z)-meson 
semileptonic decays have been made by CLEO-c at the -0(3770) 0, H, 0, II]- Absolute ex- 
clusive semileptonic branching fractions were measured for nine final states that included a 
single K, tt, K*, p, or u meson. Summing the exclusive semileptonic branching fractions gives 
EB(D° xcl ) = [6.3±0.2(stat)±0.2(syst)]% andJ2B(D+ cl ) = [15.2 ± 0.3(stat) ± 0.4(syst)]%. 
CLEO has also measured the absolute branching fractions of inclusive semileptonic de- 
cays and found B(D° -> Xe + v e ) = (6.46 ± 0.17 ± 0.13)% and B(D + -> Xe + v e ) = 
(16.13 ± 0.20 ± 0.33)% jfjj]. These measurements are consistent with, but larger than, the 
sum of the exclusive semileptonic branching fractions. Although the possibility of additional 
semileptonic modes of the D° and D + with large branching fractions is excluded, a window 
for additional semileptonic decay modes remains. 

Semileptonic decays of a D meson with an rj or 7/ in the final state have not yet been 
observed. Their discovery would open future opportunities to obtain information about 
77 — v( mixing 0]. They also probe the composition of the 77 and 7/ wave functions when 
combined with measurements of the corresponding D s semileptonic decays j§| and gauge the 
possible role of weak annihilation in the corresponding D s -meson semileptonic decays j8j. 
The process D + —>■ <pe + v e is not expected to occur in the absence of mixing between the u 
and (f> @. 

We report herein the first observation and absolute branching fraction measurement of 
D + — > r]e + v e , and results of searches for D + — ► r]'e + v e and D + — ► <pe + v e . (Throughout this 
Letter charge-conjugate modes are implied.) The data sample used for these measurements 
consists of an integrated luminosity of 281 pb _1 at the t/>(3770) resonance, and includes 
about 8 x 10 5 D + D~ events [l(|. The data were produced in e + e~ collisions at the Cornell 
Electron Storage Ring (CESR-c) and collected with the CLEO-c detector. This is the same 
data set used in Refs. [3, H, @, [Tot ]. 

CLEO-c is a general-purpose solenoidal detector. The charged particle tracking system 
covers a solid angle of 93% of 47r and consists of a small-radius six-layer low mass stereo wire 
drift chamber concentric with, and surrounded by, a 47-layer cylindrical drift chamber. The 
chambers operate in a 1.0 T magnetic field and achieve a momentum resolution of ~0.6% at 
p =1 GeV/c. The main drift chamber provides specific-ionization (dE/dx) measurements 
that discriminate between charged pions and kaons. Additional hadron identification is 
provided by a Ring-Imaging Cherenkov (RICH) detector covering approximately 80% of 
47r. Identification of positrons and detection of neutral pions and eta mesons relies on an 
electromagnetic calorimeter consisting of 7800 cesium iodide crystals and covering about 
93% of 47T. The calorimeter achieves a photon energy resolution of 2.2% at =1 GeV and 
5% at 100 MeV. The CLEO-c detector is described in detail elsewhere llll . 



The technique for our analysis was first applied by the Mark III collaboration [12| at 
SPEAR. The presence of two D ± mesons in a D + D~ event allows a tag sample to be defined 
in which a D~ is reconstructed in a hadronic decay mode. A sub-sample is then defined in 
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which a positron and a set of hadrons, as a signature of a semileptonic decay, are required 
in addition to the tag. The tag yield can be expressed as N tag = 2A r £ ) £>£>tagetag, where Ndd 
is the produced number of D + D~ pairs, £> tag is the branching fraction of hadronic modes 
used in the tag sample, and e tag is the tag efficiency The yield of tags with a semileptonic 
decay can be expressed as iV tagiSL = 2N ,DD$tag^SL e tag,SL where £>sl is the semileptonic decay 
branching fraction, including subsidiary branching fractions, and e tag ,sL is the efficiency of 
finding the tag and the semileptonic decay in the same event. From the expressions for jV tag 
and iVtag,sL we obtain 

„ N tagt st e tag -ZVtag.SL/e , , 

SL ~ ~n — ; — ~ ~iv — ' 

iv tag Ctag.SL iv tag 

where e = e tag) sL/etag is the effective signal efficiency. The branching fraction determined 
by tagging is an absolute measurement. It is independent of the integrated luminosity and 
number of D + mesons in the data sample. Due to the large solid angle acceptance and high 
segmentation of the CLEO-c detector and the low multiplicity of the events e ta g,sL ~ etagCsLj 
where esL is the semileptonic decay efficiency. Hence the ratio etag.sLAtag is insensitive to 
most systematic effects associated with the tag mode and the absolute branching fraction 
determined with this procedure is nearly independent of the tag mode. 

Candidate events are selected by reconstructing a D~ tag in one of the following six 
hadronic final states: Kgn~ , K + tt~tt~, K^tt'tt , K + tt~tt~tt° , Kgir~TT~TT + , and K + K~tt~. 
Tagged events are selected based on two variables: AE = £/£> — -Ebeain> the difference between 
the energy of the D~ tag candidate (Ed) and the beam energy (i^beam), and the beam- 



constrained mass Mb c = yE^^/c 4 — \pd\ 2 /c 2 , where po is the measured momentum of the 
D~ candidate. Selection criteria for tracks, tt° and K% candidates used in the reconstruction 



of tags are described in Ref. [L3J. If multiple candidates are present in the same tag mode, 
one candidate per tag charge with the smallest AE is chosen. The yield of each tag mode, 
and the combined yield of all six tag modes, are obtained from fits to the Mb c distributions, 
where the signal shape includes the effects of beam energy smearing, initial state radiation, 
the line shape of the ■0(3770'), and reconstruction resolution, and the background is described 
by an ARGUS function [l4j ]. which models combinatorial contributions. The data sample 
comprises approximately 163,000 reconstructed charged tags (Tabled]). 

TABLE I: Tag yields of the six D~ hadronic modes with statistical uncertainties. 



Tag Mode iV, 



D- 






11469 


± 116 


D- 


-> K+7T- 


'IT 


79933 


± 293 


D- 






25243 


± 212 


D- 


-> K+7T- 




23733 


± 185 


D- 




7T~7T + 


16446 


± 177 


D- 


-> K~K 


+ 7T- 


6785 


± 97 



All Tags 163057 ± 483 



After a tag is identified, we search for a positron and a set of hadrons recoiling against 
the tag. (Muons are not used because the CLEO-c muon identification system has poor 
acceptance in the momentum range characteristic of semileptonic D decays at the ^(3770).) 
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Positron candidates are selected based on a likelihood ratio constructed from three inputs: 
the ratio of the energy deposited in the calorimeter to the measured momentum (E/p), 
dE/dx, and RICH information [3|. Furthermore, candidates must have momenta of at least 
200 MeV/c and satisfy | cos^l < 0.90, where 9 is the angle between the positron direction 
and the beam axis. The efficiency for positron identification has been measured primarily 
with radiative Bhabha events. In the kinematic region used in this analysis, the efficiency 
rises from about 50% at 200 MeV/c to 95% just above 300 MeV/c and is roughly constant 
thereafter. The rates for misidentifying charged pions and kaons as positrons, averaged over 
the momentum range, is approximately 0.1%. Bremsstrahlung photons are recovered by 
adding showers within 5° of the positron that are not matched to other particles. 

Hadronic tracks produced in semileptonic D + decay are required to have momenta above 
50 MeV/c and | cos#| < 0.93. Pion and kaon candidates are required to have dE/dx mea- 
surements within three standard deviations (3a) of the expected value. For tracks with 
momenta greater than 700 MeV/c, RICH information, if available, is combined with dE/dx. 
The efficiencies (95% or higher) and misidentification rates (a few per cent) are determined 
with charged pion and kaon samples from hadronic D decays. 

We select 7r° candidates from pairs of photons, each having an energy of at least 30 MeV, 
and a shower shape consistent with that expected for a photon. A kinematic fit is performed 
constraining the invariant mass of the photon pair to the known 7r° mass. The candidate 
is accepted if the unconstrained invariant mass is within 3a of the nominal 7r° mass, where 
o (typically 6 MeV/c 2 ) is determined for that candidate from the kinematic fit, and the 
kinematic parameters for the 7r° determined with the fit are used in further reconstruction. 
We reconstruct 77 candidates in two decay modes. For the decay 77 — > 77, candidates are 
formed using the same procedure as for 7r° except that a ~ 12 MeV/c 2 . For 77 — > tt + 7t~tt° we 
require that the invariant mass of the three pions be within 12 MeV/c 2 of the known 77 mass. 
We reconstruct 7/ in the decay mode 7/ — > -n + n~ri. We require [m^^-^ —m n >\ < 10 MeV/c 2 . 
We reconstruct <fi candidates in the mode (f> — > K + K~ requiring \rriKK — m^l < 13.5 MeV/c 2 . 
For both the rj' and (j) these mass cuts correspond to ±3<r. 

The D~ tag and D + semileptonic decay are combined if they account for all tracks in 
the event. Semileptonic decays are identified with U = E miss — c|p miss |, where -E m i SS and 
Pmiss are the missing energy and momentum of the D + meson. If the decay products have 
been correctly identified, U is expected to be zero, since only a neutrino is undetected. 
The resolution in U is improved by constraining the magnitude and direction of the D + 



momentum to be po+ = \J i? 2 eam /c 2 — c 2 m 2 Dl and pr>+ = —pr>- [3j, respectively. Due to the 
finite resolution of the detector, the distribution in U is approximately Gaussian, centered 
at U = with a ~ 10 MeV (the width varies by mode). The number of events with multiple 
candidates varies by mode, ranging from zero for D + — > 0e + z/ e to 83% for D + — > 7/(7^77, 77 — > 
7T7T7T°)e + i> e . To remove multiple candidates in each semileptonic mode one combination is 
chosen per tag mode per tag charge, based on the proximity of the invariant masses of the 
7T°, 77, 7/, or candidates to their expected masses. 

The U distributions for D + — > f]e + u e for each 77 decay mode with all tag modes com- 
bined are shown in Fig. [TJ The U distributions in data for D + — > r]'(7i + 7T^7])e + i , e and 
D + — ► <j)(K + K~)e + v e have no entries in —250 < U < 250 MeV and are not included in Fig. 



1 [l5j. The yield for D + — > r]e + v e is determined from a binned likelihood fit to the U distri- 



bution where the signal is described by a modified Crystal Ball function with two power-law 
tails [l6| which account for initial- and final-state radiation (FSR) and mismeasured tracks. 
The signal parameters are fixed with a GEANT-based Monte Carlo (MC) simulation 17J] in 



5 



0121007-001 




> 

LU 



0.2 -0.1 0.1 0.2 




E . -dp . l(GeV) 

miss nniss v ' 



FIG. 1: Fits to the U distributions in data (filled circles with error bars) for D + — > r]e + v e for two 
r\ decay modes: (a) r/ — > 77, (b) 77 — > 7r + 7r~7r°. The solid line represents the fit of the sum of the 
signal function and background function to the data. The dashed line indicates the background 
contribution. 



fits to the data. The background functions are determined by a simulation that incorporates 
all available data on D meson decays. The backgrounds are small and arise mostly from mis- 
reconstructed semileptonic decays with correctly reconstructed tags. For D + — > rj{^)e + v e) 
there is background from D + — > Ks(7T°7r°)e + i> e and D + — > 7r°e + z/ e . The background shape 
parameters are fixed, while the background normalizations are allowed to float in fits to the 
data. The signal yields A^ tagiSL are given in Table [TTJ The fits describe the data well. After 
increasing the backgrounds by one statistical a, the probabilities that background fluctua- 
tions account for the signals are 2 x 10~ 10 and 6 x 10~ 8 , corresponding to 6.3a and 5.3a, for 
D + — > r/(77)e + z/ e and D + — > ?7(7r + 7r~7r )e + z/ e , respectively. This is the first observation of 
D + — ► 7]e + v e . 

The absolute branching fractions and 90% confidence level (C.L.) upper limits in Table HT1 
are obtained from Eq. ([T]). The signal efficiencies, e, are determined by MC simulation, 
and have been weighted by the tag yields shown in Table [H Due to differences between 
the simulation and data, corrections are applied to A^ tagiSL for 7r° and 77 finding, positron 
identification, and charged 7r and K identification. The corrections are given in Table IIHI 
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We consider the following sources of systematic uncertainty and give our estimates of their 
magnitudes in Table IIHI Where these estimates vary by decay mode a range is given. The 
uncertainty in the track finding efficiency is estimated with missing mass techniques applied 
to the data and simulation [10]. The correction to the tt° — > 77 (rj — > 77) detection efficiency 



and its uncertainty are estimated with the same method. The momentum- dependent cor 
rection to the positron identification, and its uncertainty, are obtained from comparisons of 
the detector response to positrons from radiative Bhabhas in data and MC simulations 0]. 
The corrections to the charged pion and kaon identification efficiencies, and the associated 
uncertainties, are estimated using hadronic .D-meson decays. The uncertainty in the number 
of tags is estimated by using alternative signal functions in the fits, and by counting. The 
uncertainty in modelling the background shapes in the fits to the U distributions has contri- 
butions from the uncertainties in the simulation of the positron and hadron fake rates. The 
uncertainty associated with the requirement that there be no additional tracks in tagged 
semileptonic events is estimated by comparing fully reconstructed DD events in data and 
simulation. The uncertainty associated with the shape of the signal function is estimated 
by using alternative signal functions. The uncertainty in the semileptonic reconstruction 
efficiencies due to imperfect knowledge of the semileptonic form factors is estimated by 
varying the form factors in the corresponding well-studied transitions D — > K/tt/ K*e + u e , 
by their uncertainties and assuming that the resulting changes in efficiency are a measure 
of the uncertainty in efficiencies for the semileptonic transitions being studied here 0. The 
uncertainty associated with the simulation of FSR and bremsstrahlung in the detector ma- 
terial is estimated by varying the amount of FSR modelled by the PHOTOS algorithm 18 
and by repeating the analysis with and without recovery of FSR photons. The uncertainty 
associated with the simulation of initial-state radiation (e + e~ — > DD^y) is negligible. These 
estimates are added in quadrature to obtain the total systematic uncertainties reported in 
Table [IT] 

For the upper limits in Table [Til systematic uncertainties are incorporated by combining 
with statistical uncertainties in quadrature and increasing the upper limit on the number 
of observed events by one a of the combined uncertainty. Systematic uncertainties for all 
modes are much smaller than statistical uncertainties. 

Our branching fractions for D + — ► r]e + u e and upper limits for D + — ► r]'e + v e and D + — > 
<pe + v e are compared to previous upper limits 0] in Table. [UJ The branching fractions 
measured using the two rj decay modes are consistent, and the combined branching fraction 
is consistent with the PDG upper limit. Theoretical predictions for B(D + — > 7]e + v e ) and 
B(D + — > r]'e + v e ) in the ISGW2 model 0] and a model (FK) which combines heavy-quark 
symmetry and properties of the chiral Lagrangian [l9j|, are also listed in Table. (TTJ Our 
B(D + — > r]e + i' e ) is consistent with these predictions, although our statistical uncertainty 
is large. The upper limits for D + — > r]'e + v e and D + — > <fie + v e are about two orders of 
magnitude more restrictive than previous experimental limits. 

In summary, we have made the first observation of D + — > r]e + i , e and measured the 
branching fraction, which is found to be consistent with model predictions. The sum of 
exclusive semileptonic branching fractions of the D + is approximately 0.9% smaller than 
the measured inclusive semileptonic branching fraction. Our B(D + — > r]e + v e ) increases 
the exclusive sum by about 0.1%, indicating that further exclusive rare semileptonic modes 
may await discovery. We have searched for the decays D + — ► 7]'e + v e and D + — ► <pe + v e 
and set si gnifi cantly improved upper limits for each mode. The predictions for B(D + — > 
rj'e + v e ) 0, Jil] are similar in magnitude to our upper limit. If these models are correct, an 
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observation of D + — > f]'e + u e is likely in the near future. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the A. P. Sloan Foundation, 
the National Science Foundation, the U.S. Department of Energy, and the Natural Sciences 
and Engineering Research Council of Canada. 
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TABLE II: Signal efficiencies, yields (90% C.L. intervals), and branching fractions (90 % C.L. upper limits) for D + — » r/e + ^ e , (D + — > r( e r v e 
and 0e + z^ e ) in units of xl0~ 4 in this work and comparisons to PDG [a] and two theoretical predictions ISGWII 0] and FK 19]. In the 
fourth column the first uncertainty is statistical and the second systematic. In other columns the uncertainty is statistical or statistical and 
systematic combined in quadrature. The efficiencies include subsidiary branching fractions from PDG 



Decay Mode 


e(%) 


^tag, SL 




^sl(PDG) 


B SL (ISGWII) 


#sl(FK) 




15.13 ±0.07 


32.6 ±6.7 


13.2 ±2.3 ±0.6 








D + — > r)(7r + ir~ir°)e + v e 


5.99 ±0.04 


13.3 ±4.0 


13.6 ±3.7 ±0.5 








D + — > r]e + v e (Combined) 






13.3 ±2.0 ±0.6 


<70 


11 


10 


D + — > r/' '(ir + tt~ r] — 1 7r + 7r~77)e + i^ e 


3.38 ±0.02 


(0.00, 2.30) 


<4.4 








D + -> ?/(vr + vr"r/ -> 2(7r + 7r^)7r°)e + z/ e 


0.85 ±0.01 


(0.00, 2.30) 


<17.3 








D + — > rj'e^Uf, (Combined) 


4.23 ±0.02 


(0.00, 2.30) 


<3.5 


<110 


5 


1.6 


L>+ -» <j)(K + K-)e + v e 


8.97 ±0.10 


(0.00, 2.30) 


<1.6 


<209 







TABLE III: Summary of corrections and systematic uncertainties to Bsl in percent (%). A minus 
sign indicates the MC simulation has higher efficiency than the data. 



Source Correction (%) Uncertainty (%) 



1 racking 






U.o - 1.0 


V — > 77 reconstruction 


— 


.0 


4.U 


7T nilQlIlg 


— O.o - 


— o.u 


z.o - O.O 


Electron identification 


-1.2 - 


-1.0 


1.0 


Background shape 






0.0 - 1.8 


Hadron identification 


-6.9 - 


-2.2 


0.2 - 0.8 


Number of D tags 






0.5 


Veto of unused tracks 






0.3 


Signal Shape 






0.0 - 0.4 


Simulation of form factors 






1.0 - 3.0 


Simulation of FSR 






0.6 


MC statistics 






0.9 - 1.8 


Total 


-12.6 - 


-7.6 


3.8 - 4.8 
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